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ABSTRACT 


(Distribution  Limitation  Statement  B) 

A  similitude  analysis  is  conducted  and  procedures  are 

presented  for  using  model  tests  to  determine  the  extent  of  bomb 
crater  damage  to  rigid  runways.  The  radius,  depth,  and  volume 
of  true  craters,  as  well  as  the  extent  of  concrete  cracking,  can 
all  be  determined;  however,  different  methods  to  scale  the  energy 
release  are  required,  dependent  upon  a  crater  or  camouflet  being 
formed.  Experimental  test  data  using  energy  releases  ranging 
from  5  grams  to  589  pounds  (750-pound  bomb)  of  C-4  are  used  to 
demonstrate  the  validity  of  this  model  analysis.  In  addition,  em¬ 
pirical  equations  have  been  curve-fitted  to  the  pi  terms  and  can  be 
used  for  predicting  mode  of  response,  true  crater  size,  and  extent 

of  concrete  cracking. 
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SECTION  I 


INTRODUCTION 

Whenever  a  bomb  penetrates  a  runway  or  a  charge  is  deto¬ 
nated  beneath  a  paved  surface,  the  resulting  explosion  creates  a 
crater  and  cracks  the  pavement  at  large  distances  beyond  the  crater 
The  crater  may  be  a  large  cavity  which  is  open  to  the  atmosphere, 
in  which  case  we  say  the  response  has  been  in  the  cratering  mode, 
or  it  may  remain  covered,  in  which  case  v.e  say  the  response  has 
been  in  the  camouflet  mode.  Irrespective  of  the  mode  of  response, 
rapid  repairs  require  a  knowledge  of  crater  size  and  extent  of 
concrete  cracking. 

Various  researchers  have  concerned  themselves  with 
the  somewhat  analogous  problem  of  predicting  apparent  crater  size 
from  buried  explosive  sources.  Unlike  these  earlier  studies,  con¬ 
cerned  largely  with  explosive  excavation,  we  are  not  interested 
in  apparent  crater  size.  Tree  crater  size  is  a  much  more  mean¬ 
ingful  measure  of  damage  to  us.  The  debris  which  accumulates 
in  the  crater  of  a  bomb-damaged  runway  may  have  to  be  removed, 
at  least  compacted,  before  repairs  can  commence.  In  addition, 

if  one  wishes  to  repair  runways,  the  extent  of  concrete  cracking 

,  ,  -> 

is  of  extreme  interest.  Although  procedures  exist  for  estimating 

r 

apparent  crater  size  essentially  no  accurate  procedures  have  been 
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developed  for  estimating  true  crater  size.  The  initial  intent  of 
this  study  was  to  review  the  bomb  cratering  data  that  presently 
exist  and  determine  if  model  tests  could  be  conducted  as  a  means 
for  obtaining  additional  data  inexpensively.  We  not  only  determined 
how  models  should  be  tested  but  also  devised  empirical  methods 

for  estimating  true  crater  size  and  the  extent  of  concrete  cracking. 

* 

Test  data  in  a  University  of  New  Mexico  report  '  on  charges 
ranging  in  size  from  5.  0  grams  to  589  pounds  of  C-4  appear  to 
confirm  the  applicability  of  any  empirical  equations. 

Section  II  of  this  report  deals  with  nomenclature.  What 
is  meant  by  modes  of  failure  and  more  exact  definitions  of  response 
parameters  are  presented  in  greater  detail  so  that  all  readers 
and  the  author  have  common  concepts  whenever  various  terms  are 
used.  To  evaluate  this  model  analysis,  experimental  true  crater 
data  and  pavement  cracking  data  were  required.  Section  III  de¬ 
scribes  all  test  sites  used  by  others  (we  performed  no  experimental 
tests  ourselves)  and  indicates  how  any  scaled  graph  should  read 
so  the  test  site,  the  pavement  thickness,  and  the  weight  of  charge 
is  known  at  a  glance  for  any  data  point.  Section  IV  presents  a 
detailed  model  analysis.  Assumptions  are  enumerated  and  their 

-  -v ,  *'  .  .  .  v  - 

implications  discussed  in  detail. 


*  Superscript  numbers  in  parentkeses  refer  to  references  at  tke 
end  of  this  report. 


An  evaluation  of  this  model  analysis  is  presented  in  Sec¬ 
tion  V.  Section  VI  contains  a  description  of  how  model  tests  should 
be  scaled  in  the  cratering  and  in  the  camcuflet  modes.  Under 
certain  conditions  which  are  illustrated  in  Section  VI,  models  can¬ 
not  be  used  to  simulate  prototype  response  because  the  mode  of 
response  will  change  from  cratering  in  the  prototype  to  camouflet 
in  the  model.  In  Section  VII  a  list  of  five  recommendations  for 
future  experiments  is  presented.  We  recommend  obtaining  addi¬ 
tional  data  to  establish  the  line  separating  the  cratering  and  camou¬ 
flet  modes  over  a  wider  range  in  initial  conditions,  collecting 
additional  data  over  a  wider  range  of  charge  weights  in  the  camou¬ 
flet  mode,  performing  some  model  experiments  to  study  the  in¬ 
fluence  of  concrete  strength  on  surface  related  phenomena  in 
both  response  modes,  measuring  the  seismic  velocity  of  soil  at 
the  various  test  sites,  and  perhaps  performing  some  model  tests 
of  cratering  from  bomb  detonations  inside  the  pavement.  A  sum¬ 
mary  concludes  this  report  and  is  presented  in  Section  VIII. 


SECTION  II 


NOMENCLATURE 


Before  beginning  this  model  analysis,  or  any  analysis, 
the  meanings  of  certain  expressions  must  be  clarified,  Figure  1 
and  this  discussion  are  presented  to  illustrate  what  we  mean  when 
certain  terms  are  used.  Whenever  an  explosive  charge  is  deto¬ 
nated  beneath  the  ground, a  cavity  or  void  is  formed  within  the  soil. 
If  the  energy  release  is  relatively  close  to  the  surface,  the  cavity 
or  void  vents  to  the  atmosphere  and  a  crater  is  formed.  Large 
amounts  of  ejecta  are  flung  upwards  and  outwards  if  the  cavity 
vents  to  the  atmosphere.  Some  of  this  ejecta  falls  back  into  the 
cavity,  whereas  other  amounts  of  it  settle  on  the  lip  of  the  crater. 
Because  large  quantities  of  ejecta  settle  into  the  cavity,  two  dif¬ 
ferent  craters  can  be  discussed;  the  apparent  crater  formed  by 
the  surface  of  the  ejecta  and  the  true  crater  formed  by  the  crater 
boundaries  without  regard  to  the  ejecta.  Illustrations  shown  graphi¬ 
cally  ir.  Figure  1  depict  the  differences  in  true  and  apparent  craters 
In  this  study,  we  are  concerned  with  true  crater  dimensions  only 
since  those  who  repair  bombed  runways  may  have  to  remove  the 
very  loosely  compacted  ejecta  before  filling  any  void.  In  Refer¬ 
ence  1  (which  is  the  source  for  much  of  the  experimental  data  used 
in  this  study),  Kvammen,  Pichumani  and  Dick  treat  craters  which 
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Figure  1.  Modes  of  Response  snd  Nomenclature 
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vent  or  are  open  to  the  atmosphere  as  two  separate  and  distinct 
modes  of  response  dependent  upon  whether  the  true  craters  are  some¬ 
what  hemispherical  or  cylindrical  in  shape.  Although  in  Figure  1 
we  show  croc  sections  of  both  types  of  craters  discussed  by  Kvam  - 
nen,  Pichumani  and  Dick,  we  make  no  distinction  between  ex¬ 
perimental  results  yielding  either  hemispherical  or  cylindrical 
true  craters.  In  this  report,  any  void  venting  to  the  atmosphere 
as  in  the  top  two  illustrations  in  Figure  1  is  said  to  have  responded 
in  the  cratering  mode. 

If,  on  the  other  hand,  the  energy  release  is  relatively  deep 
within  the  ground,  the  cavity  or  void  is  created  without  any  appreci¬ 
able  venting  to  the  atmosphere.  This  mode  of  response  as  shown 
in  Figure  1  is  termed  the  camouflet  mode.  Very  little  ejecta  is 
thrown  into  the  air  in  the  camouflet  mode.  A  large  spherical  void 
is  formed  beneath  the  surface  of  the  soil.  Because  physical  phe¬ 
nomena,  which  are  insignificant  in  determining  true  crater  size 
in  the  cratering  mode,  become  important  in  the  camouflet  mode, 
we  study  cratering  and  camouflet  response  as  separate  problems. 

The  responses  of  interest  in  either  the  cratering  or  camou¬ 
flet  modes  are  associated  with  true  crater  dimensions  and  the 
extent  of  circumferential  cracking  in  the  concrete.  Long  radial 
cracks  extend  way  beyond  the  circumferential  cracks  in  a  damaged 
concrete  pavement;  however,  it  is  the  area  of  concrete  covered 
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by  both  radial  and  circumferential  cracks  that  require  immediate 
repair.  Kvammen,  Pichumani,  and  Dick  in  Reference  1  present 
bomb  cratering  damage  to  runway  experimental  data  through  the 
use  of  six  parameters.  Because  we  rely  so  heavily  on  their  com¬ 
pilation  of  data  in  this  study,  we  use  these  six  measures  of  respons 
as  defined  in  Reference  1  and  as  described  below. 

(1)  True  Crater  Depth  (D)  -  Depth  from  the  bottom 

of  the  true  crater  to  the 
top  of  the  original  con¬ 
crete  surface. 

(2)  Surface  Crater  Radius  (R^)  -  The  average  radius 

of  the  crater  at  the 
top  of  the  original 
concrete  surface. 

This  average  is  cal¬ 
culated  by  taking 
the  square  root  of 
the  area  at  the  crater 
surface  divided  by 
pi. 

The  concrete  cracks  beyond 
the  surface  crater  radius. 

An  average  radius  of  circum¬ 
ferential  cracking  was  cal¬ 
culated  by  taking  the  square 
root  of  the  area  enclosed  by 
circumferential  cracking  in 
concrete  divided  by  pi. 


(4)  .Maximum  Crater 

Radius  (R^)  -  Maximum  radius  of  true  crater 

at  any  depth.  For  the  cratering 
mode  of  response  R^  usually 


(3)  Radius  of  Concrete 
Cracking  (R^) 


equals  R^.  This  parameter 
has  more  meaning  in  the  camou- 
flet  mode  of  response  where  the 
maximum  radius  occurs  far  be¬ 
low  the  surface  of  the  soil. 

(5)  True  Soil  Crater 

Volume  (v^  )  -  Volume  of  the  soil  crater  from 

the  bottom  of  the  concrete  slab 
to  the  true  crater  walls. 

(6)  Soil  and  Concrete 

Crater  Volume  (V^)  -  Volume  of  concrete  ejecta 

plus  V  .  The  volume  of 
the  crater  from  the  top  of 
the  slab  to  the  true  crater 
walls. 

Other  parameters  used  in  this  study  are  inputs  rather  than 

responses.  They  include  the  density  of  soil  p  ,  the  density  of 

6 

concrete  n  ,  the  ultimate  strength  of  the  concrete  p  ,  ,  the 
c  ult 

seismic  wave  velocity  through  the  soil  c,  the  acceleration  of  grav¬ 
ity  g,  the  depth  of  burial  for  the  explosive  charge  d,  and  the  equiv¬ 
alent  energy  release  W.  Depth  of  burial  is  always  measured 
from  the  top  of  the  original  concrete  surface  f  j  the  c.  g.  of  the 
explosive  source.  Because  different  chemical  explosives  can  be 
used  as  an  energy  source,  we  always  consider  the  energy  release  W 
as  an  equivalent  weight  relative  to  an  explosive  called  C-4.  We 

use  the  energy  of  explosion  as  our  criterion  for  relating  explosive 

(  2) 

sources.  The  energy  of  explosion  for  C-4  equals  1165  calories/ 
gram. 


SECTION  III 


•SOURCES  FOR  TEST  DATA 

The  data  used  for  all  evaluations  of  bomb  crater  damage  to 
runways  came  from  a  single  source  by  Kvammen,  Pichumani,  and 
Dick^  ^  .  This  report  contains  a  compilation  of  experimental  test 
data  from  pavement  cratering  studies.  Primarily  all  of  the  data 
in  their  report  came  from  several  different  series  of  experiments 
at  one  of  three  different  locations.  Various  experiments  at  these 
three  different  locations  furnish  us  with  enough  data  for  evaluating 
the  effects  of  charge  weight,  soil  type,  and  pavement  thickness  on 
bomb  crater  damage  to  runways. 

The  Fort  Sumner  test  site,  located  in  Southeastern  New 
Mexico,  consisted  of  a  75-feet  wide  by  1,800-feet  long  by  7-inches 
thick  concrete  slab  overlying  a  silty  sand  medium.  Blowcount,  dry 
density,  moisture  content,  specific  gravity,  and  the  Atterberg  limits 
were  obtained  at  five  different  test  boring  locations  and  three 
different  trench  sites.  Table  I  from  Reference  l  gives  an  indica- 
tion  of  soil  conditions  at  Fort  Sumner  by  showing  the  log  for  one 
of  the  bore  holes.  Additional  sor.  details  should  be  obtained  from 
Appendix  I  in  Reference  1.  No  base  or  subbase  underlies  the  7- 
inches  of  paving  at  Fort  Sumner.  The  concrete  was  poured  directly 
onto  the  soil;  a  runway  paving  technique  typical  of  old  World  War  II 


Table  I 


FORT  SUMNER  SOIL  PROPERTIES 


Depth, 

ft 

Unified 

Classification 

Description 

Vd» 

pcf 

g 

LL 

PI 

w, 

% 

Station  SI  6 

2 - 

3 - 

SC 

Clayey  Sand 

21 

9 

_ 

16.0 

4 - 

n 

1 

■ 

SM 

Silty  Sand 

104. 62 

2.67 

16 

2 

8.  5 

o 

8 - 

SM 

Silty  Sand 

104.67 

2.  69 

I 

9.5 

1 

1 

1  1 

SP-SM 

Sand  w/Silt 

121. 14 

2.66 

E 

D 

1  1 

■ 

■ 

l  c - 

SW -SM 

Gravelly  Sand 
w/Silt 

110.27 

2.66 

E 

a 

H 

1  3 

1 

1 4 

CH 

Fat  Clays 

114.34 

2.67 

59 

36 

24.0 

_ 

* 
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vintage  airfields.  Because  of  this  paving  technique,  the  thickness 
of  the  concrete  was  approximately  7  inches  with  significant  varia¬ 
tions  in  thickness  occuring.  Most  of  the  concrete  pads  were  20- 
feet  long  by  10-feet  wide;  however  a  few  were  20-feet  long  by  15- 
feet  wide.  The  average  compressive  strength  of  the  concrete  at 
Fort  Sumner  was  10,400  psi.  Additional  concrete  test  data  for 
Fort  Sumner  may  also  be  obtained  from  Appendix  I  in  Reference  1. 

The  Hays,  Kansas,  test  site  consisted  of  a  150-feet  wide 
by  5,  600-feet  long  strip  of  concrete  overlying  a  clay  soil.  Con¬ 
crete  thickness  ranged  from  11 -inches  for  the  center  pads  to  8- 
inches  for  the  outlying  pads.  With  the  exception  of  blow  count , 
the  same  joil  properties  were  measured  at  Hays  as  at  Fort  Sumner. 
One  of  the  boring  logs,  as  presented  in  Reference  1,  is  shown  in 
Table  II  to  give  an  indication  of  soil  conditions  at  Hays,  Kansas. 

A  wider  variation  in  soil  properties  does  exist  at  Hays  than  at 
Fort  Sumner.  No  gravel  base  course  underlies  the  paving  at 
Hays  either.  Each  concrete  pad  at  Hays  was  approximately  20-feet 
long  by  12.  5-feet  wide.  The  average  compressive  strength  of  the 
concrete  at  Hays  was  10, 100  psi,  almost  identical  to  the  compres¬ 
sive  strength  for  the  concrete  at  Fort  Sumner.  Additional  soil 
and  concrete  information  about  the  Hays  test  site  may  be  obtained 
from  Appendix  II  in  Reference  1 . 


Table  11 


Depth, 

ft 

Unified 

Classification 

Description 

Yd« 

pcf 

Dari  Brown  Clay  95.60  2.68  71  45  28.0 


CH 

Gray  Brown  Clay 

102.  60 

2.  6 

Light  Brown  Clay  97.70  2.70  60  37  28.0 


Light  Silty  Clay  108.90  2.  67  32  13  1  5.5 


Light  Silty  Clay  115.80  2.67  43  26  16.5 

w /Large  Propor¬ 
tion!  of  Caliche 


In  addition  to  the  Fort  Sumner  and  Hays  cratering  experi¬ 
ments,  the  University  of  New  Mexico  generated  extra  cratering 
data  in  a  series  of  what  they  termed  prototype  and  model  tests^. 
Only  data  from  two  prototype  tests  are  included  in  Reference  1. 
These  tests  were  conducted  on  a  14-inch  thick  concrete  slab  which 
has  been  spread  over  a  6-inch  thick  gravel  base.  Under  the  con¬ 
crete  and  gravel  was  a  6-foot  thick  clay  subgrade.  The  prototype 
pavement  was  divided  by  joints  into  25-foot  by  25-foot  panels. 

No  information  concerning  prototype  concrete  strength  is  given 
in  Reference  1,  and  essentially  no  information  about  soil  or  gravel 
base  properties. 

The  series  of  ten  University  of  New  Mexico  model  experi¬ 
ments  is  described  in  greater  detail  than  the  prototype  experiments. 
In  the  model  tests,  a  2.75-inch  thick  concrete  slab  was  poured  over 
a  1 . 25-inch  gravel  base.  Under  the  concrete  and  gravel  lay  a 
15-ir.ch  thick  layer  of  clay.  The  model  pavement  was  divided  by 
joints  into  panels  4  feet  10  inches  square.  The  average  compres¬ 
sive  strength  of  the  model  concrete  was  3705  psi  which  is  consid¬ 
erably  lower  than  the  average  compressive  strengths  for  the  Fort 
Sumner  and  Hays  airfields.  Density  and  moisture  contents  were 
reported  for  the  model  soils;  these  averaged  112.  0  lb/ft^  and 
11. 0%,  respectively. 
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The  vast  majority  of  the  data  used  in  an  evolution  of  our 
modeling  laws  comes  from  the  Fort  Sumner  and  Hays  Test  sites. 
These  aata  are  supplemented  by  a  l<t w  University  of  New  Mexico 
prototype  and  model  experiments.  Different  sizes  of  explosive 
charge  were  used  at  various  test  sites;  5-pound  C-4,  15-pound 
C-4,  and  25-pound  C-4  charges  were  placed  at  various  depths  under 
7.  0-inch  thick  concrete  at  Fort  Sumner,  8. 0-inch  thick  concrete 
at  Hays,  and  1 1 . 0-inch  thick  concrete  at  Hays.  In  addition,  250- 
pound  (MK-81),  500-pound  (MK-82),  and  750-pound  (M-117)  bombs 
were  placed  at  various  depths  under  11.0  inch  thick  concrete  at 
Hays  only.  The  two  prototype  University  of  New  Mexico  experi¬ 
ments  used  1.  5  pound  C-4  charges  and  the  ten  model  experiments 
used  various  size  C-4  charges  ranging  from  5.  0  grams  to  17.5 
grams.  Table  III  gives  an  indication  of  charge  dimensions  and 
the  equivalent  energy  release.  This  writer  believes  that  the  three 
different  types  of  bombs  carried  Tritonal,  a  very  energic  chemi¬ 
cal  explosive  with  an  energy  of  explosion  of  1770  calories/gram^. 

Because  C-4  is  a  much  less  energic  explosive,  an  energy  of  ex- 

2 

plosion  of  1165  calories /gram  ,  the  actual  weight  of  explosive  in 
each  bomb  was  multiplied  by  1 . 52  (1770/1165)  to  obtain  an  equiva¬ 
lent  C-4  energy  release.  The  influence  of  bomb  casing  on  the 
energy  release  is  assumed  to  be  insignificant. 
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Table  III 


CHARGE  SIZES  AND  ENERGY  RELEASES 


Explosive 

Source 

Diameter  of 
Explosive 
Source(in.  ) 

Length  to 
Diameter 
Ratio 

Actual 

Weight  of 
Explosive(lb) 

Equivalent 
C-4  Energy 
Release 

5.  0  lbs  C-4 

4.  0 

1.87 

5.  0 

5 

15.0  lbs  C-4 

6.  0 

1. 66 

15.  0 

1  5 

25.  0  lbs  C-4 

7.0 

1.75 

25,  0 

25 

250  lb  Bomb 

9.0 

5.48 

101. 5 

154 

500  lb  Bomb 

10.  75 

6. 

192.  5 

292 

750  lb  Bomb 

16.  1 

3.  1  ! 

387.  0 

589 

Notice  also  that  the  bombs  are  many  diameters  longer  than  the 
C-4  plosive  charges.  We  are  of  the  opinion  that  this  length  to 
diameter  aspect  ratio  has  an  insignificant  influence  on  true  crater 
dimensions  and  the  extent  of  concrete  cracking;  nevertheless,  we 
mention  this  observation  for  the  sake  of  completeness.  All  depths 
of  burial  were  measured  and  will  be  used  in  our  analysis  relative 
to  the  c.g.  of  the  explosive  source.  No  geometric  charge  dimensions 
were  given  for  the  University  of  New  Mexico  prototype  and  model 
explosive  sources;  hence,  this  information  cannot  be  included  in 
Table  IU.  The  prototype  and  model  experiments  did  report  their 
charge  weights,  and  the  explosive  used  in  these  experiments  was 
C-4. 
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Throughout  the  many  graphical  comparisons  which  are  about 
to  be  made  for  evaluations  of  modeling  laws,  we  use  experimental 
test  data  from  the  various  test  sites  and  numerous  different  sixes 
of  explosive  charges.  To  tell  for  any  given  data  point  what  was 
the  charge  size,  what  was  the  thickness  of  the  overlying  concrete, 
and  how  large  was  the  energy  release,  we  had  to  devise  a  syste¬ 
matic  procedure  for  presenting  the  test  data  from  Reference  1. 

Table  IV  illustrates  this  procedure  graphically.  Basically  the 
procedure  involves  use  of  shading  to  denote  the  test  site  and  con¬ 
crete  thickness,  and  use  of  symbol  shape  to  denote  magnitude 
of  the  energy  release. 


The  white  symbols  are  for  data  points  taken  at  Fort  Sumner,  the 
black  symbols  are  for  data  points  taken  at  Hays  with  11.0  inches  of 
concrete  as  overburden,  and  the  speckled  or  gray  figures  are  for 
data  points  taken  at  Hays  with  8.  0  inches  of  concrete  as  overbur¬ 
den.  All  model  tests  irrespective  of  the  size  of  the  energy  release 
are  shown  by  the  symbol  X,  and  the  two  prototype  1.  5-pound  ex¬ 
periments  are  denoted  by  a  plus  sign  +.  The  larger  the  energy 
release,  the  more  sides  are  given  to  a  symbol  (a  circle  is  assumed 
to  have  an  infinite  number  of  sides).  Table  IV  shows  that  the  energy 
release  increases  as  the  symbols  change  from  triangles,  to  diamonds 
to  squares,  to  pentagons,  to  hexagons,  and  finally  to  circles. 

All  graphs  in  this  report  for  evaluating  model  approaches  to  ex¬ 
perimental  bomb  crater  damage  to  runways  use  the  system  of 
symbols  presented  in  Table  IV. 

Most  experiments  report  true  crater  radius  at  the  surface 
R  ,  average  radius  of  concrete  cracking  F  ,  maximum  true  crater 

L  M 

radius  at  any  depth  R^,  depth  of  true  crater  D,  volume  of  true 
crater  without  including  concrete  V^ ,  and  volume  of  true  crater 
and  damaged  concrete  V^.  In  a  few  experiments,  some  of  the  para¬ 
meters  are  not  reported  so  we  cannot  include  them.  For  example, 
all  of  the  bomb  data  report  surface  effects  R  and  R,;  however, 
we  cannot  include  in  this  st  :dy  D,  V^,  V^  and  for  some  of  the 
bomb  experiments  because  some  craters  were  never  excavated. 
Reference  1  also  does  not  report  radii  R^  and  R ,  in  the  prototype 


University  of  New  Mexico  tests.  We  use  all  data  whenever  it  is 
reported.  In  making  evaluations  using  the  University  of  New 
Mexico  prototype  and  model  test  data,  we  assume  that  the  thick¬ 
ness  of  the  pavemen*-  is  that  of  the  concrete  slab,  and  that  the 
gravel  base  is  part  ot  the  soil  media. 

In  some  of  our  graphical  comparisons,  we  use  cratering 
data  for  experiments  with  no  paving  for  confinement.  These  data 
come  from  a  compendium  compiled  at  Waterways  Experiment 
Station  '  ,  and  are  used  to  demonstrate  that  under  certain  condi¬ 

tions  the  scaled  thickness  of  pavement  and  scaled  concrete  strength 
are  insignificant  parameters  in  determining  certain  true  crater 
dimensions.  We  do  not  use  different  symbols  or  shadings  to  denote 
test  site  location  or  magnitude  of  the  energy  release  for  data  from 
Reference  3,  cratering  without  confinement  from  paving;  neverthe¬ 
less,  these  true  cratering  data  come  from  tests  in  clay,  sand, 
moist  loess,  and  wet  silt,  and  include  equivalent  energy  releases 
of  0.25,  0.50,  1.0,  8.0,  25,  27,  54,  64,  256,  320,  3560,  40,000, 
and  320,  000  pounds  of  TNT.  Tables  1  and  2  in  the  Appendices  of 
Reference  3,  and  in  Reference  3  itself,  should  be  studied  by  those 
seeking  additional  information  about  these  data. 
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SECTION  IV 


MODEL  ANALYSIS 

Any  model  analysis  reflects  the  definition  of  the  problem. 
Physical  phenomena  that  are  not  considered  cannot  appear  in  the 
resulting  pi  terms.  In  this  study  we  exercise  our  opinion  and 
engineering  judgment  to  list  parameters  which  we  feel  may  be  sig¬ 
nificant.  A  problem  such  as  this  one  to  define  true  crater  dimen¬ 
sions  and  predict  radius  of  concrete  cracking  is  very  difficult  and 
subject  to  differences  in  professional  opinion.  In  defining  this 
problem,  we  attempt  to  indicate  what  assumptions  are  being  made. 
Subsequent  nondimensional  plots  are  presented  to  illustrate  the 
degree  to  which  this  model  analysis  is  valid  and  to  provide  physi¬ 
cal  insight  into  the  significance  of  the  analysis. 

The  parameters  which  we  consider  include  the  total  energy 
release  in  the  buried  explosive  charge  W  and  the  depth  of  burial 
from  the  top  surface  of  the  concrete  d.  We  do  not  include  para¬ 
meters  defining  the  geometry  of  the  explosive  source,  nor  do  we 
list  parameters  for  studying  the  explosive  process  and  the  propaga¬ 
tion  of  shocks  within  the  explosive.  Two  main  reasons  exist  for 
deleting  these  details  from  the  analysis.  First,  we  wish  to  keep 
our  list  of  parameters  as  small  as  is  practically  possible.  Second, 
we  are  not  interested  in  the  intensity  of  pressures  and  other  lo- 
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calized  details  in  the  vicinity  of  the  charge;  we  are  interested  in 
gross  overall  response  on  a  much  larger  scale,  such  as  true  crater 
dimensions.  By  listing  only  W  and  d  as  input  parameters  related 
to  the  explosive  source,  we  are  assuming  that  the  energy  release 
is  instantaneous  relative  to  cratering  times  and  originates  from 
a  source  whose  size  is  insignificant. 

The  major  difficulty  in  this  problem  is  the  selection  of 
appropriate  parameters  for  characterizing  soil  behavior.  Out  of 
necessity  we  assume  that  the  soil  is  homogeneous  to  an  infinite  depth, 
isotropic,  and  a  single  phase  medium.  Admittedly  pore  air  and 
pore  water  pressures  within  the  granular  structure  of  the  soil  in¬ 
fluence  shear  strength;  however,  this  process  is  so  poorly  under¬ 
stood  in  soil  dynamics  that  we  cannot  expect  to  reflect  complex 
behavior  in  an  applied  problem  seek  ng  numerical  values.  We  char¬ 
acterize  the  properties  of  soil  using  three  parameters;  the  density 
of  the  media  P  ,  the  seismic  velocity  of  the  soil  c  and  the  ac¬ 
celeration  of  gravity  g.  These  three  parameters  are  selected  be¬ 
cause  this  investigator  has  had  some  success  in  applying  them 

(4) 

to  a  related  problem  ,  the  apparent  crater  dimensions  caused 
by  the  detonation  of  a  buried  explosive  charge.  Selecting  only  three 
parameters  for  characterizing  soil  does  permit  the  ground  to  behave 
as  a  fluid  close  to  the  charge  and  as  a  solid  at  greater  standoff 
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distance.  Inertial  effects  are  incorporated  into  the  analysis  through 

the  density  P  •  Gravitational  effects  are  included  because  the  ac- 
’  s 

celeration  of  gravity  g  is  listed,  as  well  as  P  .  The  seismic 

s 

velocity  parameter  c  when  considered  in  combination  with 

2 

represents  the  soil's  constitutive  strength.  The  product  P  c 

s 

is  a  measure  of  the  compressibility  of  the  ground.  It  takes  on 
fim.e  numerical  values  for  all  soil  types  and  never  equals  zero. 
Other  parameters  (unconfined  compressive  strength,  cone  indexes, 
etc.  )  which  might  be  selected  as  a  measure  of  constitutive  strength 
often  are  numerically  either  infinite  or  zero  for  certain  extremes 
in  soil  conditions.  These  parameters  are  unacceptable  measures 
of  constitutive  strength  as  they  lead  to  indeterminate  solutions 
even  when  we  are  physically  aware  that  craters  do  form.  By  se¬ 
lecting  only  three  soil  parameters,  we  represent  the  maximum 
number  of  three  physical  phenomena  (gravitational,  inertial,  and 
constitutive  effects)  with  the  absolute  minimum  number  of  para¬ 
meters  . 

We  do  not  include  strain  rate  effects  in  this  analysis  by 
listing  a  rate  of  strain  parameter  because  we  believe  these  effects 
are  insignificant.  Certainly  all  cratering  events  are  at  high  rates 
of  strain,  but  the  parameters  P#  and  c  are  measures  of  constitu¬ 
tive  behavior  at  quite  high  rates  of  strain  also.  Provided  con¬ 
stitutive  strength  is  not  appreciably  modified  by  changes  in  the 
rates  of  strain  between  the  smallest  and  the  largest  charges,  this 


assumption  is  acceptable.  Usually  many  orders  of  magnitude  change 
in  rate  of  strain  are  required  to  appreciably  modify  constitutive 
properties.  Going  in  charge  weight  from  1  pound  to  1000  pounds 
while  holding  everything  else  constant  approximately  changes  the 
rate  of  strain  by  only  one  order  of  magnitude,  an  insignificant  amount 
for  appreciably  modifying  constitutive  properties. 

Because  the  ground  is  covered  by  slabs  of  concrete,  we 
must  list  parameters  for  defining  concrete  properties.  We  assume 
that  the  concrete  has  little  or  no  reinforcing  steel,  and  that  it 
extends  to  infinity  in  the  plane  of  the  earth's  surface.  Although 
we  list  no  steel  reinforcing  properties  in  this  analysis,  we  present 
our  views  about  simulating  pavements  with  large  amounts  of  rein¬ 
forcing  after  we  analyze  test  results  and  discuss  how  to  construct 
model  airfield  cratering  experiments.  The  concrete  pavement 
is  of  thickness  h  and  has  inertial  properties  represented  by  the 
density  parameter  and  constitutive  strength  represented  by  the 
ultimate  compressive  strength  of  the  concrete  o^.  We  assume 
strain  rate  effects  are  insignificant  in  the  concrete  as  in  the  soil, 
and  that  the  ultimate  tensile  strength  of  concrete,  toughness,  and/ 
or  other  strength  characteristics  may  be  thought  of  as  constant 
ratios  relative  to  the  ultimate  compressive  stress  0  ^  . 

So  far  in  this  analysis,  we  have  characterized  the  explosive 
source,  the  soil,  and  the  concrete  pavement.  Next  we  must  decide 
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what  responses  are  of  interest.  Each  of  these  responses  or  results 
is  a  function  of  the  input  explosive,  soil,  and  concrete  parameters; 
they  are  not  functions  of  one  another.  For  any  given  experimental 
test,  we  would  like  to  know  whether  the  response  will  be  in  the 
cratering  or  the  camouflet  modes.  Mode  shape  is  a  nondimensional 
parameter  which  we  shall  call  "MODE1'.  We  would  like  to  predict 
true  crater  size,  what  quantity  of  concrete  becomes  ejecta,  and 
over  how  large  a  region  the  concrete  is  cracked.  Because  Kvammen, 
Pichumani,  and  Dick  in  Reference  1  report  true  crater  depth  D, 
surface  crater  R^,  average  radius  of  circumferential  concrete 
cracking  R^,  maximum  crater  radius  R^,  true  soil  crater  volume 
Vj,  and  soil  and  concrete  crater  volumes  VJt  we  use  these  six 
convenient  dependent  parameters  as  our  response  parameters. 

From  these  six  parameters,  one  can  estimate  true  crater  size 
and  the  extent  of  damage  to  the  pavement  which  requires  repair. 

Table  V  lists  the  eight  independent  input  parameters,  plus 
the  seven  dependent  response  parameters  which  we  have  just  de¬ 
scribed,  and  presents  their  fundamental  dimensions  in  the  engi¬ 
neering  system  of  force  F,  length  L,  and  time  T. 
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Table  V 

LIST  OF  PARAMETERS 

Independent  Input  Parameters 

Fundamental 

Symbol 

Pa  rameter 

Dimensions 

W 

Energy  Release  in  Explosive 

FL 

d 

Depth  of  Explosive  Burial 

L 

P 

s 

Density  of  Soil 

ft2/l4 

g 

Acceleration  of  Gravity 

L/T2 

c 

Soil  Seismic  Velocity 

L/T 

p 

c 

Density  of  Concrete 

ft2/l4 

a 

ult 

Ultimate  Strength  Concrete 

F/L2 

h 

Thickness  of  Concrete 

L 

Dependent  Response  Parameters 

Symbol  Parameter 

MODE  Mode  of  Response 

D  Depth  of  Crater 


Fundamental 

Dimensions 


L 


l 

I 


f 


R^  Surface  Crater  Radius  L 

R  Radius  of  Concrete  Cracking  L 

Ct 


Because  there  are  fifteen  parameters  and  three  fundamental 
dimensions,  we  can  develop  five  independent  pi  terms  and  seven 
dependent  response  pi  terms  by  applying  the  Buckingham  Pi  Theo¬ 
rem.  The  mathematics  of  obtaining  these  pi  terms  are  straight¬ 
forward  and  will  not  be  repeated  in  this  report;  no  new  assump¬ 
tions  are  introduced  into  this  analysis  because  we  go  through  the 
tedious  algebra  of  obtaining  these  nondimensional  ratios.  If  one 
solves  for  pi  terms  with  respect  to  the  parameters  P  ,  c,  and 
d,  he  will  obtain  the  pi  terms  shown  in  Table  VI. 


Table  VI 


PI  TERMS  -  OBTAINED  BY  SOLVING  RELATIVE  TO  P,  c,  and  d 

s 


Independent  Pi  Terms 


Dependent  Response  Pi  Terms 


TTj  =  h/d 

n  _  =  P  /p 
Z  c  s 

n  /  2 

n  ,  =o  ../  P  c 
3  ult  s 

v  4-  gd/c2 

TT  =  W /  P  C2  d3 
5  s 


"  ,  =  MODE 
o 

"  7-Vd 

n  8  =  R2/d 

"  9  =  Vd 

n  l  o  s  D/d 

"u  1  V"3 

n12  =  V2/d3 


f 
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Pi  terms  may  be  multiplied  or  divided  by  one  another 
to  form  new  pi  terms.  They  may  be  taken  to  any  power.  These 
manipulations  are  performed  for  the  convenience  of  the  investi¬ 
gator.  Although  the  pi  terms  in  Table  VI  are  a  perfectly  acceptable 
set,  this  writer  finds  it  convenient  to  rewrite  pi  term  4  as  the 
ratio  of  energy  release  relative  to  the  energy  expended  in  over¬ 
coming  gravitational  effects  rather  than  use  it  in  its  present  for¬ 
mat  of  gravitational  effects  relative  to  constitutive  effects.  Di¬ 
viding  pi  term  5  by  pi  term  4  to  obtain  a  new  pi  term  4,  then  taking 
the  1/4  root  of  the  new  pi  term  4,  1/3  root  of  pi  term  5,  1/3  root 
of  pi  term  11,  and  1/3  root  of  pi  term  12  yields  the  pi  terms 
shown  in  Table  VII.  The  pi  terms  shown  in  Table  VII  are  used 
in  all  subsequent  analyses. 


Table  VII 

PI  TERMS  -  USED  IN  ALL  EVALUATIONS 


Independent  Pi  Terms 

Dependent  Response  Pi  Terms 

tt  =  h/d 

tt  =  MODE 

0 

TT  =  o  /  o 

TT  =  R  /d 

2  c  s 

7  1 

TT  =  a  /  P  C2 

Tt  =  R  /d 

3  ult  8 

8  2 

TT  =  W 1/4  /  P  1/4  gl/4d 

"  =  R  ,/d 

4  S 

9  3 

TT  *  Wl/3/  P  1/3  C2/3  d 

5  s 

"io‘D/d 

"u-v  ,l/3'd 

-  V  ,/3/d 
"  12  •  V2  /d 
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Provided  that  the  parameters  which  we  have  listed  as 
measures  of  concrete  strength  and  soil  strength  are  appropriate 
and  provided  that  constitutive,  inertial,  and  gravitational  effects 
are  the  physical  phenomena  involved  in  the  cratering  event,  Equa¬ 
tion  (1)  for  any  of  the  response  pi  terms  in  Table  VII  is  a  predic¬ 
tion  equation  interrelating  all  variables  in  a  bomb  crater  damage 
to  runway  analysis. 


n  dependent 
parameter 


i  (  h 
'T* 


ult 


W 


1/4 


W 


1/3 


(1) 


D  c 

8 


TiTjrr 


This  equation  defines  a  six-parameter  space  of  nondimensional 

products.  It  is  much  too  large  a  space  to  empirically  generate 

an  entire  solution  for  Equation  (1)  as  presented.  Fortunately,  as 

a  first-order  approximation,  certain  parameters  may  be  treated 

as  constant,  and,  as  a  result,  Equation  (1)  can  be  simplified. 

We  will  treat  0  0,  g,  c,  and  n  w  as  constants  in 

c  s  ult 

subsequent  analysis.  The  fact  that  the  first  three  of  these  five 

parameters  are  essentially  constant  should  disturb  no  one.  Very 

little  variation  of  lees  than  4  15%  occurs  in  P  ,  p  ,  and  g  from 

c  s 

experiment  to  experiment.  The  other  two  parameters,  c  and 
,  do  have  the  capability  for  significant  variation.  Unfortu- 
nately  seismic  velocity  c  is  not  reported  in  the  test  data  compiled 
in  Reference  1.  For  most  sands  and  clays,  c  will  range  from 
around  600  to  6000  fpe.  We  will  assume  c  is  constant  for  the  soils 
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used  in  Reference  1  because  we  cannot  assign  numerical  values 

(4) 

until  this  parameter  is  measured.  Earlier  work  on  an  analo¬ 
gous  problem,  apparent  crater  size  from  buried  explosive  charges, 
indicates  that  treating  c  as  a  constant  for  various  types  of  soil 
does  not  significantly  modify  apparent  crater  geometry. 

For  all  practical  purposes,  0  was  a  constant  for  most 
of  the  tests  reported  in  Reference  1.  At  the  Fort  Sumner  test  site, 
the  average  ultimate  compressive  strength  from  twelve  tests  on 
concrete  cylinders  equaled  10,400  psi,  and,  at  Hays,  this  same 
test  on  nine  concrete  cylinders  gave  0  as  10,  100  psi.  The 
vast  majority  of  all  data  was  taken  at  Fort  Sumner  and  Hays. 

In  the  University  of  New  Mexico  model  experiments,  the  average 
concrete  ultimate  strength  was  3705  psi,  considerably  lower  than 
at  Hays  and  Fort  Sumner.  The  ultimate  concrete  strength  is  not 
reported  in  the  two  University  of  New  Mexico  prototype  tests. 

We  wilt  assume  0  ^  was  a  constant  in  all  experimental  tests,  but, 
whenever  model  test  results  disagree  with  Fort  Sumner  and  Hays 
results,  one  should  consider  if  the  cause  of  any  disagreement  could 

be  low  values  of  0  ,  . 

ult 

If  we  assume  p  ,  p  ,  g,  c,  and  a  are  constants,  we 
c  s  ult 

can  treat  these  parameters  as  abstract  numbers,  and  write  a 
simplified,  four -parameter  version  of  Equation  (1).  Equation  (2) 


2,3 


-LuU'L-.-aWSP 


is  this  simplified  version  of  Equation  (1).  Note  that  r>  /  P  and 

c  s 

2 

n  /pc  disappear  from  Equation  (1),  because  these  terms 
ult  s 

are  now  invariant,  and  that  two  of  the  independent  quotients  in 

Equation  (2)  are  now  dimensional  because  portions  of  those  pi 

1/4 

terms  have  been  treated  as  abstract  numbers.  Because  W  /d 
1/3 

and  W  /d  are  now  dimensional,  we  use  W  in  equivalent  pounds 
of  C-4  and  d  in  feet  of  burial. 

,h  w1/4  wu\ 

n  dependent  ‘  1  >  ~3T  1  (Z' 

parameters 

Equation  (2)  with  a  few  variations  depending  upon  mode  of 
response  and  response  parameter  being  studied  is  used  for  evalu 
ating  this  model  analysis. 


SECTION  V 


EVALUATION  OF  MODEL  ANALYSIS 

Mode  of  Response 

Equation  (2)  indicates  that  the  mode  of  response  (cratering 
mode  or  camouflet  mode)  should  be  a  function  of  normalized  con¬ 
crete  thickness  ~  and  two  different  normalized  charge  weights 

W 1  /3  wl/^4 

— - —  and  — - —  .  Actually  mode  sh«_pe  is  very  insensitive  to 
d  d 

normalized  pavement  thickness,  so  Equation  (3)  can  be  used  in¬ 
stead  of  Equation  (2)  for  predicting  mode  of  response. 


MODE  =  f 


W 
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) 


(3) 


mode  y  d  d 

To  prove  that  Equation  (3)  is  appropriate,  consider 
the  experimental  data  presented  in  Figure  2.  All  data  points  shown 
in  Figure  2  follow  our  standard  convention  for  plotting  data. 

In  addition,  either  a  vertical  slash  or  a  horizontal  slash  bisects 
each  data  point  in  Figure  2.  If  the  slash  is  vertical  on  this  plot 


W 


1/3 


of  ?. .  0  +  log  — ~ — r  versus  2.0+  log 


W 


1/4 


,  the  response  was 


in  the  cratering  mode,  and  if  the  slash  is  horizontal,  the  response 
was  in  the  camouflet  mode.  The  number  two  was  arbitrarily  added 
to  each  data  point  because  the  writer  did  not  care  to  have  negative 
values  for  the  logarithm  of  any  normalized  parameters.  This 
manipulation  only  translates  the  X  and  Y  axis  in  any  plot;  it  does 
not  modify  conclusions.  All  logarithms  in  this  figure  and  subse- 
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Response 


quent  figures  are  to  the  base  ten.  The  charge  weight  W  in  Equa¬ 
tion  (3)  and  Figure  2  is  in  equivalent  pounds  of  C-4,  while  the  depth 
of  burial  d  is  in  feet  of  burial. 

Notice  in  Figure  2  that  a  curved  line  separates  the  data 
responding  in  the  cratering  mode  from  the  data  responding  in  the 
carnouflet  mode.  If  the  data  fall  to  the  left  and  above  this  line, 
the  response  was  in  the  carnouflet  mode,  and,  if  the  data  fall  be¬ 
low  and  to  the  right  of  this  line,  the  response  was  in  the  cratering 
mode.  A  few  data  points  have  both  a  solid  slash  and  a  dashed 
slash;  the  mode  for  the  response  of  these  points  is  uncertain, 
but  the  mode  represented  by  the  solid  slash  is  the  more  probable 
of  the  two  modes.  Essentially  all  data  fail  on  the  appropriate  side 
of  the  line  separating  the  two  modes;  only  a  few  exceptions  exist. 
These  exceptions  are  close  to  the  line,  and  one  must  remember 
that  we  are  treating  P  ,  g,  and  especially  c  as  abstract  numbers 

3 

in  the  parameters  represented  by  the  X  and  Y  coordinates  in 
Figure  2.  The  parameter  ~  for  the  data  in  Figure  2  ranges  from 
a  low  of  0.  06  to  a  high  of  1 . 39.  A  value  of  —  greater  than  1 . 0 
means  that  the  c.g.  of  the  charge  was  actually  up  in  the  concrete 
pavement.  Hence,  the  data  in  Figure  2  range  through  a  few  charges 
which  were  actually  up  in  the  pavement,  to  many  data  points  for 
charges  buried  over  sixteen  pavement  thicknesses.  Further  response 
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evaluations  will  reveal  that  responses  in  the  cratering  mode  are 


independent  of  pavement  thickness  and  that  all  responses  e:  cept 
perhaps  the  surface  effects  (R^/d  and  R  /d)  are  independent  of 
pavement  thickness  in  the  camouflet  mode.  If  other  responses 
are  independent  of  pavement  thickne  ss,  then  intuitively  one  should 
feel  that  mode  of  response  would  ct  independent  of  pavement  thick¬ 
ness. 

The  curve  which  we  passed  between  the  data  points  in 

Figure  2  has  an  equation  given  by: 

W1^4  5  W1^3 

(2.0+  log  — - —  )  =  1 .  3981  (2.0  +  log  — - — )  tanh  (2.0+  log  — - — ) 

Equation  (4)  is  a  functional  method  of  expressing  Equation  (3). 

1/3 

If  after  substituting  for  W  / d  in  Equation  (4)  the  calculated 

1/4 

value  of  W  /d  is  greater  than  the  intended  experimental  value, 

then  the  response  falls  in  the  camouflet  mpde.  Otherwise  if  the 

1/4 

calculated  value  of  W  /d  is  less  than  the  intended  experimental 

value,  the  response  should  fall  in  the  cratering  mode.  Equation 

1/3 

(4)  is  only  valid  for  values  of  2. 0  +  log  W  / d  from  1.50  to  1.85 
because  this  equation  is  a  curve-fit  to  a  limited  quantity  of  experi¬ 
mental  data,  all  of  which  fall  within  a  narrow  band  on  a  plot,  as 
shown  in  Figure  2. 

Either  Figure  2  or  Equation  (4)  may  be  used  to  predict 
mode  of  response.  Mode  of  response  is  important  because  the 


(4) 


crater  dimensions  and  extent  of  concrete  cracking  scale  in  differ¬ 
ent  manners  dependent  upon  whether  the  response  is  in  the  camou- 
flet  or  cratering  modes. 

Cratering  Mode 

In  the  cratering  mode,  all  true  crater  dimensions  are  in¬ 
dependent  of  all  independent  normalized  parameters  in  either 

1/3 

Equations  (1)  or  (2),  except  W  /d.  This  observation  means  that 
the  shape  of  the  true  crater  depends  upon  neither  gravitational 
effects  nor  the  strength  and  inertial  effects  of  the  concrete  pave¬ 
ments.  In  other  words,  the  explosive  energy  release  is  so  over¬ 
powering  that  the  pavement  provides  insignificant  confinement, 
and  gravitational  effects  only  influence  the  free  fall  of  the  ejecta, 
an  event  which  could  not  possibly  modify  true  crater  dimensions. 

To  demonstrate  that  true  crater  dimensions  are  a  function 
1/3 

only  of  W  / d  in  the  cratering  mode,  Figures  3  through  6  present, 

respectively,  scaled  plots  of  surface  crater  radius,  crater  depth, 

soil  crater  volume,  and  soil  plus  concrete  crater  volumes  as  func- 
1/3 

tions  of  W  /d.  The  data  uaed  in  these  figures  come  for  all  the 
Reference  1  teat  sites,  and  the  explosive  charges  range  from  a 
few  grams  up  through  bombs.  The  parameter  h/d  varies  from 
a  low  of  0.  07  to  a  high  of  1. 39.  Seven  data  points  represent  test 
conditions  with  the  c.g.  of  the  explosive  source  up  in  the  concrete 
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Figure  3.  Normalized  Surface  Crater  Radius  in  Cratering  Mode 


Figure  5.  Normalized  Soil  Crater  Volume  in  Cratering  Mode 


Normalized  Soil  Plus  Concrete  Crater  Volumes  in  Cratering  Mode 


(1.00£h/d<  1.39).  Figures  3  through  6  are  log  plots,  with  the 
number  2.  0  arbitrarily  added  to  the  numerical  value  of  a  test  result, 
so  negative  values  would  not  occur  when  the  logarithm  was  less 
than  one.  As  in  Figure  2,  this  manipulation  shifts  the  X  and  Y 
axes,  but  does  not  modify  conclusions.  Some  scatter  occurs  in 
the  data  presented  in  Figures  3  through  6;  however,  this  observer 
feels  that  it  is  random  rather  than  systematic.  All  logarithms  have 
been  taken  to  the  base  ten  in  these  true  crater  dimension  figures; 
therefore,  a  difference  of  two  in  test  results  is  represented  by 
a  factor  of  approximately  0.  3  in  Figures  3  through  6  (logarithm  of 

2.  0  equals  0.  30103).  The  range  for  all  the  data  in  these  figures  is 
significantly  less  than  a  factor  2.  0  in  the  response  pi  term.  Figure 

3,  Rj/d,  has  the  largest  scatter  among  any  of  the  true  crater 
figures,  but  this  parameter  should  scatter  more  because  relative 
to  the  other  cratering  response  parameters  R^/d  is  poorly  defined. 
Often  the  slope  of  the  true  crater  near  the  surface  asymptotically 
approaches  the  horizon;  hence,  different  observers  easily  assign 
various  numerical  values  to  R^  when  they  cannot  tell  exactly 
where  a  crater  begins.  The  parameter  R^  has  not  been  presented 
as  a  measure  of  crater  size  because  in  the  vast  majority  of  cases 

equals  R^  for  responses  in  the  cratering  mode. 

Straight  lines  have  been  least-square  curve-fitted  to  the 
data  in  Figures  3  through  6.  The  equations  to  these  lines  which 


functionally  relate  the  response  pi  terms  to 


/d  are  given  by: 


R 

d 

D 

d 

V 

V 


-  =  2.  1  55  /W 
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0.865 


=  2. 312 


W 


1/3) 


0. 6833 


1/3  ; 

=  2.  046  * 


1/3 


0. 7849 


1/3  / 

V-  =  2.372  * 
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0.8227 


(5) 

(6) 

(7) 

(8) 


In  addition,  standard  deviations  have  been  calculated  by  dividing 

theoretically  predicted  responses  from  Equations  (5)  through  (8) 

into  experimentally  observed  responses.  This  manipulation  which 

yields  new  distributions  about  a  mean  o.f  unity  permits  the  compu- 

1/3 

tation  of  standard  deviations  for  R^/d,  D/d,  /d,  and 

1/3 

V^  /d  about  their  respective  predictive  equation.  Table  VIII 

summarizes  these  results.  As  indicated  in  Table  VIII,  one 

standard  deviation  essentially  equals  12%  except  for  the  less  well 

conditioned  a^/d  where  it  equals  20.4%.  This  observer  believes 

these  results  are  excellent.  Further  improvements  in  accuracy 

2 

might  be  achieved  if  one  had  actual  values  of  P  c  to  substitute 

8 

into  the  independent  pi  term;  however,  practical  problems  would 
require  some  averaging  procedure  for  obtaining  an  effective  value 
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of  P  c  because  of  changes  in  0  c  with  depth.  In  the  opinion 

8  S 

of  this  writer,  which  presently  cannot  be  substantiated  with  experi- 

2 

mental  data,  variations  of  c  with  depth  in  real  soil  will  pre¬ 
vent  orders  of  magnitude  improvements  in  accuracy. 


Table  VIII 

STANDARD  DEVIATION  CRATER  SIZE  IN 
CRATERING  MODE 


Parameter 


R1  /d 
D/d 


Standard  Deviation 
20.4% 

12.2% 

12.7% 

U.9% 


As  additional  evidence  that  h/d  is  an  insignificant  para¬ 
meter  in  the  cratering  mode,  plot*  have  been  made  from  data  in 
Reference  2  of  true  crater  surface  radius,  true  crater  depth, 
and  true  crater  volume  without  any  pavement  present  to  act  as 
constraint.  Figures  7,  8  and  9  which  present  these  plots  com¬ 
pare  this  unconstrained  true  crater  data  to  the  least-square  curve- 
fit  lines  and  the  standard  deviations  about  these  lines  from,  respec¬ 
tively,  Figure  3  for  surface  crater  radius.  Figure  4  for  crater 
depth,  and  Figure  6  for  soil  and  concrete  volume.  The  airfield 
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Figure  7.  Normalized  True  Crater  Radius 


'Vwiiin.'nn  '<1  hw  *sd>t 


Figure  8.  Normalized  True  Crater  Depth  Without  Pavement 


Normalized  True  Crater  Volume  Without  Pavement 


cratering  experiments  appear  to  predict  cratering  results  for 
h/d  equal  to  zero,  especially  when  one  considers  the  scatter  in 
both  runway  cratering  results  as  represented  by  the  3cr  standard 
deviation  bounds  anu  in  unconstrained  ciaters  as  depicted  by  in¬ 
dividual  data  points.  Without  the  presence  of  pavement,  surface 
crater  radius  is  slightly  larger,  but  crater  depth  is  slightly  shal¬ 
lower  especially  for  shallow  depths  of  burial  than  would  be  predicted 
with  the  presence  of  pavement,  Crater  volumes  are  essentially 
identical  with  or  without  the  presence  of  pavement.  These  obser¬ 
vations  substantiate  further  that  true  crater  dimensions  are  influ¬ 
enced  only  slightly  in  the  cratering  mode  by  normalized  pavement 
thickness.  If  pavement  thickness  does  not  influence  crater  size, 
this  observer  does  not  believe  that  pavement  strength,  the  pi  term 

0  i  /  will  have  any  influence  on  true  crater  geometry. 

2 

The  parameter  a  P  c  equaled  zero  when  no  pavement  was 
present  and  equaled  some  finite  value  in  Reference  1  experiments  ; 
hence,  true  crater  size  should  a1  so  be  independent  of  this  scaled 
parameter  in  the  cratering  mode. 

So  far,  we  have  evaluated  cratering  parameters;  we  have 
not  determined  average  radiuw  of  circumferential  cracking  in  the 
concrete,  a  surface  phenomenon  which  is  influenced  by  normalized 


concrete  strength,  o  /  P^c  .  Figure  10  presents  R^d  as  a 
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KANSAS  I  KANSAS  I  NEW  MEXICO 


Normalized  Radius  of  Concrete  Cracking  in  Cratering  Mode 


function  of  W  /d  for  all  cratering  mode  data  in  Reference  1. 

This  figure  indicates  that  normalized  radius  of  concrete  cracking 

1/3  2 

is  a  function  of  W  /d,  is  a  function  of  a  /  p  c  ,  but  is  not 

VI  It  s 

a  function  of  ~  .  Most  of  the  data  in  Figure  10  collapse  into  a 
d 

1/3 

unique  function  relating  R^/d  *°  ^  / d.  Many  different  values 

of  h/d  ranging  from  0.  07  to  1.  39  are  included  in  these  data,  even 
though  separate  plots  are  not  being  constructed  for  different  values 
of  h/d.  Therefore,  normalized  radius  of  cracking  is  independent 
of  h/d  within  the  limits  being  studied  in  h/d.  The  three  data  points 
which  depart  the  most  from  test  results  in  Figure  10  are  all  for 
University  of  New  Mexico  model  tests.  In  the  model  experiments, 
the  average  ultimate  concrete  strength  equaled  3705  psi;  whereas, 
in  all  other  experiments,  the  average  ultimate  concrete  strength 
was  within  a  few  hundred  psi  of  10,200  psi.  Because  the  concrete 
strength  was  weaker  in  the  model  tests  than  in  all  other  experi¬ 
ments,  the  concrete  cracked  a  greater  scaled  distance  in  model 
experiments.  Obviously  the  concrete  cracks  because  of  up-heave 
in  the  underlying  soil.  Up-heave  of  the  soil  is  not  greatly  influenced 
by  scaled  pavement  thickness  or  scaled  pavement  strength;  however, 
a  scaled  amount  of  up-heave  will  crack  a  weak  concrete,  whereas 
a  stronger  concrete  may  resist  cracking  in  the  same  scaled  amount 
of  up-hea\e. 

A  least-squares  fit  to  the  data  in  Figure  10  may  be  achieved 
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using  a  straight  line.  This  equation  is  only  valid  for  concretes 
with  an  ultimate  strength  of  approximately  10,200  psi.  The  equa¬ 
tion  is: 

0.  9a 31 

I  provided  a  =  10,200  psi  (9) 

One  standard  deviation  about  this  line  equals  16.  5%  if  the  model 
data  are  excluded  from  the  other  data.  Throughout  this  analysis, 
we  ignore  the  presence  of  any  reinforcing  and  the  size  of  slabs. 
Obviously  joints  in  stabs  can  arrest  propagating  cracks,  and  rein¬ 
forcing  can  effectively  create  a  stronger  concrete,  especially  in 
tension.  Considering  that  these  quantities  are  not  scaled  from 
test  to  test,  this  observer  is  amazed  to  find  one  standard  devia¬ 
tion  is  only  16.  5%.  One  can  scale  reinforcing  and  slab  dimensions 
if  models  are  being  made  of  a  specific  situation;  however,  handling 
these  complications  in  devising  a  formula  creates  a  much  larger 
space,  a  space  which  would  be  too  Urge  to  develop  an  entire  solu¬ 
tion  using  only  experimental  results. 


R. 


■-  4.247 


W 


1/3 


Camouflet  Mode 

In  the  camouflet  mode,  experimental  results  will  show 

that  normalized  true  crater  maximum  radius  R^/d,  true  crater 

depth  D/d,  and  soil  plus  concrete  crater  volume  V^^/d  are 

2 

all  insensitive  functions  of  h/d  and  a  /  P  c  ,  but  their  dimen- 

ult  s 
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sions  do  depend  upon  both  W^^/d  and  W*  ^/d  in  either  Equation  (1) 
or  (2).  This  observation  implies  that  the  shape  of  the  true  crater 
in  the  camouflet  mode  depends  on  both  gravitational  effects  and 
constitutive  effects  in  the  soil;  however,  the  presence  of  a  concrete 
pavement  has  only  a  minimal  influence  on  crater  size.  If  these 
ascertainments  are  correct,  one  cannot  conduct  a  mode  test  in 
the  camouflet  mode  that  rigorously  satisfies  all  pi  terms  unless 
troublesome  and  expensive  test  techniques  are  employed.  Charge 
weight  cannot  simultaneously  be  scaled  as  both  the  cube  of  the 
depth  of  burial  and  the  fourth  power  of  the  depth  of  burial.  For¬ 
tunately,  this  apparent  impass  can  be  overcome  by  employing  ap¬ 
proximate  techniques. 

(4) 

This  writer  in  earlier  studies  indicated  that  apparent 
crater  dimensions  can  be  simulated  if  the  normalized  parameters 
W^^/d  and  W^^/d  are  combined  through  an  empirical  obser¬ 
vation  to  form"\/(W*'  ^/d)*  (W^^/d)  or  W^^/d.  In  other  words, 
we  distort  both  pi  terms  by  making  one  quantity  too  large  and  the 
other  too  small  so  compensating  errors  essentially  provide  cor¬ 
rect  results.  Apparently  this  tame  empirical  observation  works 
with  reasonable  accuracy  in  scaling  true  crater  dimensions  within 
the  camouflet  mode. 

To  demonstrate  that  true  crater  dimensions  are  an  approx- 
7/24 

imate  function  of  W  / d  in  the  camouflet  mode.  Figures  1 1  through 
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13  present,  respectively,  scaled  plots  of  maximum  crater  radius 

R^/d,  crater  depth  D/d,  and  soil  plus  concrete  crater  volumes 
1  j  3  V  /  2  4 

/  /d  as  functions  of  W  /d.  The  data  used  in  these  figures 
comes  from  the  Fort  Sumner,  Hays,  and  University  of  New  Mexico 
model  experiment  test  site.  Unfortunately,  no  bomb  cratering 
experiments  were  conducted  in  the  camouflet  mode;  hence,  air¬ 
field  camouflet  cratering  experiments  only  cover  a  very  small 
range  in  size  of  explosive  charge,  from  5  to  25  pounds  of  C-4. 

The  parameter  h/d  does  not  range  through  as  wide  a  variation 
in  the  camouflet  mode  as  in  the  cratering  mode  because  the  camou¬ 
flet  mode  is  associated  with  relatively  deep  depths  of  burial. 

In  the  camouflet  mode,  h/d  appears  to  be  a  secondary  parameter 
in  determining  true  crater  size  because  the  data  in  Figures  11 
through  13  include  values  from  0.06  to  0.46  without  being  forced 
to  consider  h/d  as  significantly  modifying  results. 

As  in  the  cratering  mode,  straight  lines  have  been  least- 
squares  curve-fitted  to  the  Fort  Sumner  and  Hays  data  in  Figures  11 
through  13.  The  University  of  New  Mexico  model  experiment  data 
are  not  included  in  determining  the  equation  for  this  line,  but  these 
model  experiment  data  points  are  shown  graphically  in  Figures  11 
through  13. 
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Normalized  True  Crater  Depth  in  Camouflet  Mode 


Normalized  Soil  Plus  Concrete  Crater  Volume  in  Camouflot  Mode 


The  equation  for  these  lines  are  given  by: 


=  1. 053 
a 


r7/24 


D  ,'W7'24 

r*  2-244  l-— 


0.8548 


0. 4324 


=  1.718 


0.  8648 


Standard  deviations  were  also  calculated  by  dividing  theo¬ 
retically  predicted  responses  from  Equation?  (10)  through  (12) 
into  experimentally  observed  responses  as  seen  in  Figures  11 
through  13.  Table  IX  presents  the  standard  deviations  for  R^/d, 
D/d,  and  V^^/d  in  the  camouflet  mode  as  determined  through 
this  normalization  manipulation  for  providing  a  large  sample  with 


a  mean  of  unity. 


Table  IX 


STANDARD  DEVIATION  OF  CRATER  SIZE 
IN  CAMOUFLET  MODE 


Parameter 


Standard  Deviation 


R3/d 


13.  55, 


v2l/3/d 


1  0.  5% 

11.9% 


Ml 


One  standard  deviation  for  any  of  the  parameters  in  the  camouflet 
mode  essentially  equals  or  is  even  a  little  less  than  the  same  stan¬ 
dard  deviation  in  the  cratering  mode,  see  Table  VIII.  Figures  11 
through  13  may  appear  to  give  larger  scatter  than  their  counterparts 
in  the  cratering  modes,  Figures  3  through  5;  however,  the  scale 
to  which  they  are  plotted  is  greatly  enlarged  relative  to  the  scales 
presented  in  the  cratering  mode.  To  some  degree,  the  Fort  Sumner 
(open  symbols)  data  appear  to  fall  below  the  Hays  data  in  Figures 
3  through  5.  This  systematic  error  may  indicate  that  the  effective 
value  of  c  for  the  larger  depths  of  burial  at  Fort  Sumner  is  higher 
than  the  average  value  of  c  for  large  depths  at  Hays.  The  Hays 
and  Fort  Sumner  results  do  appear  as  parallel  lines  with  a  small 
offset.  This  observation  indicates  that  a  constant  such  as  c  and 

even  the  less  probable  P  or  g  would  improve  the  scatter  were 

8 

they  known  and  inserted  into  the  independent  pi  term. 

Additional  evidence  that  h/d  has  an  influence,  but  a  very 

small  one,  on  true  crater  size  in  the  camouflet  mode  is  presented 

in  Figures  14  and  15  where  normalized  camouflet  maximum  radius 

7/24 

and  volume  are  plotted  as  functions  of  W  / d  for  data  from 
Reference  3,  which  has  no  confining  pavement.  Figures  14  and 
15  compare  unconstrained  camouflet  data  to  the  least-squares 
curve-fitted  lines  and  their  standard  deviations  from,  respectively, 
Figure  11  for  true  maximum  camouflet  radius  and  Figure  13  for 
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Figure  15.  Normalized  Camouflet  Volume  Without  Pavement 


soil  and  concrete  volume.  Because  much  of  the  data  in  Reference  3 


are  at  even  greater  scaled  depths  of  burial  than  the  camouflet 
data  for  pavements  in  Reference  1,  the  lines  which  were  curve- 
fitted  must  be  extended.  The  dashed  lines  in  Figures  14  and  15 
are  extensions  to  the  curve-fitted  lines.  All  unconstrained  data 
lie  close  to  the  line  and  its  extension  especially  when  one  considers 
the  width  of  the  three  standard  deviation  confidence  band;  however, 
the  unconfined  data  does  fall  below  the  runway  camorf,~t  line. 

This  observation  implies  that  h/d  has  some  influence  on  true 
crater  dimensions  in  the  camouflet  mode,  but  this  influence  is 
small.  Nothing  about  similitude  theory  implies  that  the  curve  which 
we  fit  to  data  has  to  be  straight.  Because  these  additional  data 
fall  primarily  in  the  lower  left-hand  corner  of  Figures  14  and  15, 
a  curved  relationship  might  be  more  appropriate  and  represent  a 
more  accurate  curve-fit  over  a  wider  range  in  test  results.  The 
present  straight-line  relationships  should  only  be  applied  for 
W  ^  ^/d  greater  than  0.2. 

Radius  of  crater  at  the  surface  of  a  camouflet  and  radius  of 
concrete  cracking  are  phenomena  which  can  be  influenced  by  sur¬ 
face  conditions,  especially  concrete  properties.  Figure  16  is  a 
three-dimensional  graph  of  normalised  surface  crater  radius  in 

the  camouflet  mode  as  a  function  of  normalized  concrete  thickness 

7/24 

h/d  and  normalized  energy  release  W  / d.  Two-dimensional 
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graphs  of  R^/d  as  a  function  of  W  /d  exhibit  enough  extra 
scatter  to  imply  that  scaled  concrete  thickness  is  indeed  a  signifi¬ 
cant  parameter.  The  general  graphical  solution  to  a  relationship 

as  shown  in  Figure  16  would  be  a  surface  over  and  through  the 
,  7/24 

h/d  versus  W  /d  plane.  Because  most  of  the  experimental 

data  points  were  taken  within  a  narrow  band  on  the  h/d  versus 
7/24 

W  /d  plane,  Figure  16  presents  essentially  one  line  from 

among  the  infinite  number  of  lines  defining  the  surface  which  is 

a  solution  for  R^/d  as  a  function  of  h/d  and  W  /d.  Figure 

16  can  be  used  to  graphically  predict  surface  crater  radius  in  the 

7/24 

camouflet  mode  provided  h/d  and  W  /d  are  related,  as  in 

7  /24 

Figure  16.  If,  for  example,  W  /d  equaled  1.95  and  h/d 
equaled  0.75,  no  prediction  could  be  made  of  R^/d  because  (  -  per 
imental  data  have  not  been  obtained  to  define  the  surface  which 


is  the  solution  in  this  regime.  Because  Figure  16  verifies  the 

modeling,  but  does  not  present  a  complete  solution,  no  equation 

can  be  given  for  predicting  R^/d  and  no  standard  deviation  can  be 

calculated  as  an  indication  of  accuracy. 

The  other  surface  parameter  R^/d  probably  is  not  a  func- 

2 

tion  of  h/d,  but  it  may  be  a  function  of  a  Pgc  as  in  the 
cratering  mode.  Figure  17  presents  radius  of  concrete  cracking 

7 /24 

data  as  a  function  of  W  /d  for  all  values  of  h/d  (h/d  values 
in  these  experiments  actually  range  from  0.  06  to  0.  46),  Large 
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difference  in  charge  weight  are  not  given  by  the  data  in  Figure  17; 


however,  Fort  Sumner  and  Hays  data  which  have  essentially  the 

same  value  of  <7  both  appear  to  yield  the  same  solution.  The 
ult 

University  of  New  Mexico  modeling  experiments  appear  to  over¬ 
estimate  R?/d  in  the  camouflet  mode.  Remember  this  same  group 
of  experiments  also  overestimates  R^/d  *n  the  cratering  mode, 
but  0  is  significantly  lower  in  the  model  tests  than  in  the 
Fort  Sumner  and  Hays  experiments.  A  concrete  v/hich  was  weaker 
than  desired  would  result  in  greater  radii  of  concrete  cracking. 
Rather  insufficient  quantities  of  model  data  are  shown  in  Figure 

17,  only  three  data  points,  so  whether  or  not  R2^  is  also  a  func- 
2 

tion  of  a  /  0  c  cannot  conclusively  be  determined  without 
additional  tests.  An  equation  has  been  fitted  to  all  the  data  except 
the  model  data  in  Figure  17.  This  relationship  which  is  given  by 
the  following  equation  should  only  be  used  provided  that  n 
equals  approximately  10,200  psi. 


ult 


R  /  7/24  \ 

(V) 


%  1.270 


provided  o  =  10,200  psi 
ult 


One  standard  deviation  for  Equation  (13)  equals  11. 1%  so  R.,/d 
in  the  camo\flet  mode  can  be  predicted  with  as  great  if  not  greater 
accuracy  than  R^/d  1°  the  cratering  mode.  Distance  between 
joints  in  concrete  slabs  and  quantities  of  reinforcing  are  not  scaled 
for  the  data  in  Figure  17,  just  as  these  quantities  are  not  scaled 


(13) 


6) 


in  the  cratering  mode.  Apparently  the  appropriate  magnitude  of 
damage  can  be  established  even  though  these  quantities  are  dis¬ 


torted. 

1  /3 

We  have  not  studied  /d  in  the  camouflet  mode  because 

V  essentially  equals  V  in  the  camouflet  mode.  Only  a  very 
1  2 

small  additional  volume  of  concrete  is  added  to  the  soil  volume 

1/3 

in  a  camouflet.  Hence,  in  this  mode,  we  do  not  study  V  /d, 
just  as  in  the  cratering  mode  we  do  not  study  R^/d.  Neither  lends 
additional  insight  to  this  analysis. 
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SECTION  VI 


FABRICATING  MODELS 

The  principle  behind  similitude  theory  is  that  individual 
variables  do  not  have  to  be  invariant;  instead,  if  nondimensional 
ratios  of  products  and  quotients,  called  pi  terms,  remain  invari¬ 
ant,  different  systems  are  equivalent.  In  this  model  analysis, 
we  consider  the  parameters  listed  in  Table  V.  These  parameters 
do  not  he  .'*  _o  be  numerically  equal  in  a  different  system.  Provided 
that  the  independent  pi  terms  listed  in  Table  VII  are  equal  in  dif¬ 
ferent  systems,  the  dependent  response  pi  terms  listed  in  j>ble  VII 
will  be  equal.  Our  major  difficulty  in  designing  a  model  experi¬ 
ment,  which  uses  the  same  soils,  same  explosive,  and  tests  in  a 
1-g  gravitational  environment  as  exists  here  on  earth,  is  that  with 
these  constraints  pi  terms  4  and  5  from  Table  VII  cannot  be  simul¬ 
taneously  satisfied. 

If  we  are  to  use  a  physically  smaller  model  experiment 
for  predicting  true  crater  size  and  the  extent  of  cracking  in  a 
prototype  runway  system,  then  some  length  ratio  ^  will  be  assigned 
to  our  characteristic  length,  depth  of  burial  d.  Adopting  the  con¬ 
vention  that: 

.  .  _  ( _ )  model 

r  (  )  prototype 

where  any  parameter  can  be  inserted  into  the  brackets,  similitude 
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requirements  are  imposed  if  we  use  the  same  materials  in  a  smaller 
model  here  on  earth  in  a  1-g  gravitational  field.  These  requirements 
are: 

d  =  X 

r 


to  \  =1.0 

\  alT 


Pi  terms  2  and  3  in  Table  VII  are  automatically  satisfied  by  the 
proceeding  requirements.  Pi  term  1  will  be  satisfied  if  we  main¬ 
tain  geometric  similarity  and  the  h^  ratio  equals  X  .  Difficulties 

4 

arise  because  we  cannot  simultaneously  scale  W  by  X  as  re- 

3 

quired  by  pi  term  4  and  as  X  as  required  by  pi  term  5.  We  can¬ 
not  build  a  model  to  a  smaller  scale,  using  the  same  materials 
and  in  a  1-g  gravitational  field,  which  satisfies  a  completely  general 
model  as  presented  in  Equation  (1)  or  its  simplified  version. 

Equation  (2).  We  can  only  build  models  satisfying  approximate 
versions  of  Equation  (1}  or  must  seek  other  model  building  techniques. 

If  we  first  determine  the  mode  of  response,  a  model  can 
be  fabricated.  The  same  materials  could  be  used  in  smaller  model 
systems  within  the  earth's  natural  gravitational  field  provided 
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3 

charge  weight  W  is  scaled  as  X  in  the  cratering  mode  and  as 
24 /7 

X  in  the  camouflet  mode.  Scaling  charge  weights  in  these 

manners  would  satisfy  all  responses  if  our  evaluation  of  the  model 
analysis  is  correct  as  presented  in  the  preceding  section.  The 
necessity  of  scaling  concrete  strength  and  maintaining  geometric 
similarity  in  the  thickness  of  runways  could  probably  be  ignored 
if  crater  geometry  only  was  studied;  however,  being  more  precise 
and  satisfying  these  additional  conditions  should  present  no  dif¬ 
ficulties,  so  we  recommend  scaling  concrete  strength  and  main¬ 
taining  geometric  similarity  in  pavement  thickness. 

To  create  a  model  concrete,  one  might  prefer  using  pea 
gravel  or  other  small  stone  in  place  of  a  prototype  large  aggregate. 

In  addition,  the  size  of  reinforcing  rods  should  also  be  scaled  so 
that  steel  does  not  provide  a  disproportionate  amount  of  strength 
in  the  model.  These  actions  will  require  other  cement-to-water 
ratios  in  a  model  concrete  because  the  surface  areas  being  wetted 
are  no  longer  scaled  unless  the  grain  size  of  sand  and  cement 
are  also  scaled,  a  very  difficult  and  impractical  undertaking. 

Several  trial  mixes  should  lead  to  a  concrete  with  the  appropriate 
strength.  We  have  indicated  that  models  do  not  necessarily  have 
to  model  concrete  strength;  however,  we  would  point  out  that  all 
data  to  date  have  been  taken  within  a  narrow  band  on  the  mode  shape 
curve  of  W*^/d  versus  W^^/d,  Figure  2.  Experiments  conducted 
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outside  this  band  could  well  result  in  pi  terms  such  as  h/d  and 
/  2 

a  /  p  c  becoming  significant.  Any  observations  in  the  pre¬ 
ceding  section  can  only  be  said  to  rigorously  apply  for  tests  within 
the  experimental  band  shown  in  Figure  2. 

Under  some  circumstances,  model  tests  cannot  be  conducted 
in  the  cratering  mode,  although  we  believe  they  are  always  pos¬ 
sible  in  the  camouflet  mode.  As  an  example  of  why  a  model  test 
might  be  impossible  in  the  cratering  mode,  consider  the  following 
illustration.  Assume  that  one  wishes  to  build  a  1/5  scale  model 

to  determine  true  crater  size  for  a  prototype  194.4-pound  charge 

1/3 

buried  16.  3  feet  deep.  The  quantity  (2.  0  +  log  W  Id)  equals 

1/4 

1.55  and  (2.  0  +  log  W  /d)  equals  1 . 36  which  indicate  in  Figure  2 

that  the  response  is  just  within  the  cratering  mode.  If  the  response 

is  in  the  cratering  mode,  we  should  scale  W  as  ^  ,  and  our 

1/5  scale  model  yields  a  charge  weight  of  1.  555  pounds  at  a  depth 

of  3.26  feet.  Before  testing  our  model,  we  should  check  its  mode 

1/3 

of  response.  The  model  quantity  (2.  0  +  log  W  / d)  equals  1.  55 

1/4 

as  does  the  prototype,  but  the  quantity  (2.  0  +  log  W  / d)  equals 

1. 535.  Figure  2  indicates  that  the  response  of  the  model  falls 
in  the  camouflet  mode,  so  we  would  not  obtain  the  proper  results 
in  a  model  experiment.  This  example  demonstrates  that  a  model 
experiment  can  change  modes  from  the  cratering  to  the  camouflet 
mode.  We  have  not  been  able  to  devise  a  set  of  circumstances  in 
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which  a  model  would  be  in  the  cratering  mode  when  the  prototype 
was  in  the  camouflet  mode.  Generally  model  points  fall  above 
corresponding  prototype  points  and  sometimes  to  the  left  of  them 
as  shown  in  Figure  2;  hence,  it  is  difficult  to  slip  from  the  camou¬ 
flet  into  the  cratering  mode. 

One  method  does  exist  for  satisfying  all  pi  terms  in  Equation 

(l);  however,  very  special  equipment  is  required.  If  one  has  access 

to  a  centrifuge,  an  artificial  body  force  can  be  generated  which 

approximates  the  influences  of  gravity  on  a  bed  of  soil.  All  pi 

terms  in  Equation  (l)  will  be  satisfied  if  the  charge  weight  W 
3 

is  scaled  as  X  and  the  acceleration  of  gravity  is  scaled  as 
l /X  .  Because  the  factor  1 /X  will  be  greater  than  one,  this  tech¬ 
nique  will  theoretically  work.  The  weakness  in  this  approach 
is  that  the  larger  the  scale  factor,  the  longer  the  arms  on  the 
centrifuge  and  the  higher  the  angular  velocity.  In  addition,  people 
having  centrifuges  are  generally  reluctant  to  have  explosives  threat¬ 
ening  their  apparatus. 

By  way  of  a  review.  Table  X  has  been  compiled  to  summarize 
how  each  parameter  wiil  usually  be  scaled.  The  mode  of  response 
must  be  checked  for  both  model  and  prototype  conditions.  If  both 
model  and  prototype  systems  full  within  the  same  mode  of  response, 
the  scale  factors  presented  in  Table  X  should  work  fairly  well. 

Notice  only  the  charge  weight  is  scaled  differently  in  various  modes. 


Table  X 


SCALE  FACTORS  FOR  MODEL  CRATERING  EXPERIMENTS 


Parameter 

Symbol 

Scale  Factor 

In  Cratering 
Mode 

Scale  Factor 

In  Camouflet 

Mode 

Depth  Burial 

d 

X 

X 

Thickness  of  Concrete 

h 

X 

X 

Density  of  Concrete 

0 

c 

1.0 

l.  0 

Ultimate  Strength 
of  Concrete 

ault 

1 . 0 

1  .  0 

Acceleration  of  Gravity 

g 

Not  Scaled 

1 . 0 

Density  of  Soil 

p 

S- 

1  .  0 

1 . 0 

Seismic  Velocity  of  Soil 

c 

l.  0 

1  .  0 

Charge  Weight 

W 

x3 

x24/7 

Volume  of  Crater 

V 

x3 

x3 

Depth  of  Crater 

D 

X 

X 

Radius  of  Crater  or  of 
Concrete  Cracking 

R 

X 

X 

SECTION  VII 


RECOMMENDA  tions 
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in  initial  conditions.  Figure  2  and  the  equation  associated  with  it 

1/3 

should  presently  be  used  only  for  values  of  (2.  0  +  log  W  / d) 

1/3  , 

between  1 .  55  and  1 . 82,  values  o.t  W  /d  between  0.  3 55  and 

0.651.  This  line  separating  modes  of  response  goes  somewhere 

which  remains  to  be  determined  outside  this  narrow  range.  If 

a  cratering  condition  falls  in  a  region  where  the  mode  of  response 

is  not  known,  predictions  cannot  be  made.  Additional  cratering 

data  would  not  only  establish  mode  of  response  more  accurately 

but  also  provide  additional  data  on  magnitude  of  responses. 

As  a  second  priority,  we  strongly  recommend  collecting 

additional  data  in  the  camouflet  mode.  Absolutely  no  experiments 

using  bombs  were  conducted  in  this  regime,  and  the  few  model 

tests  that  were  performed  using  gram  charges  have  many  critics 

skeptical  of  the  results.  The  vast  majority  of  data  which  everyone 

appears  prepared  to  accept  involves  energy  releases  between 

5  and  25  pounds  of  C-4.  Our  conclusions  in  the  camouflet  mode 

must  be  accepted  as  tentative  because  a  factor  of  5  difference 

in  charge  weight  only  represents  a  very  6mall  60%  variation  in 

7/24 

the  important  independent  variable  W  A  few  experiments 

using  750*<pound  bombs  in  the  camouflet  mode  would  be  expensive; 

however,  they  would  increase  this  very  small  60%  variation  in 
7  /  24 

W  to  a  much  more  significant  402%  variation.  Supplement 


four  to  five  bomb  tests  with  additional  model  experiments  havipg 
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charges  less  than  a  pound  and  the  analysis  in  the  camouflet  mode 
wouid  either  be  on  more  solid  grounds  or  would  he  modified. 

As  a  third  priority,  we  recommend  model  experiments  to 
determine  the  influence  of  concrete  strength  on  surface-related 
responses  in  both  the  cratering  and  camouflet  modes.  The  vast 
majority  of  experiments  which  were  conducted  at  Fort  Sumner  and 
Hays  all  had  concretes  with  essentially  the  same  ultimate  compres¬ 
sive  strength  of  10,200  psi.  Cracking  of  concrete  in  both  modes 
and  radius  of  surface  crater  in  the  camouflet  mode  are  responses 
which  appear  to  have  a  dependency  on  concrete  strength.  Until 
additional  test  results  are  acquired,  inadequate  amounts  of  infor- 

2 

mation  exist  to  empirically  determine  how  the  pi  term  a  /  p  c 

ult  s 

influences  Rj/d  in  the  camouflet  mode  and  especially  R^/d 
in  both  the  camouflet  and  cratering  modes. 

Recommendations  of  lower  priority,  but  which  could  be 
easily  accomplished,  include  measuring  the  seismic  velocity  as 
a  function  of  depth  at  the  Fort  Sumner  and  Hays  test  sites  and 
conducting  more  experiments  to  simulate  bombs  up  in  the  con¬ 
crete.  If  the  seismic  velocity  c  were  accurately  known  as  a  func¬ 
tion  of  depth,  an  averaging  procedure  might  be  developed  for  ob¬ 
taining  an  effective  seismic  velocity  c  which  could  reduce  scat¬ 
ter.  Provided  moisture  contents  were  much  the  same  under  the 
pavements  at  Fort  Sumner  and  Hays,  measurements  of  c  taken 
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now  should  not  differ  too  greatly  from  the  actual  value  of  c  at 
the  time  of  testing. 

Only  seven  tests  were  conducted  with  the  center  of  charges 
in  the  concrete.  This  experimental  condition  does  lend  itself  to 
model  tests,  so  if  information  is  desired  cn  cratering  from  charges 
within  the  pavement  or  even  because  of  surface  burst,  we  recom¬ 
mend  using  models.  Scaled  thickness  of  concrete  h/d  and  scaled 

2 

concrete  strength  o^/  P  c  would  probably  become  significant 
parameters,  and  would  certainly  be  important  parameters  for  air 
bursts . 

If  additional  tests  as  recommended  herein  were  conducted, 
the  general  overall  need  to  obtain  some  response  data  outside  the 
band  on  the  W^^/d  versus  W^^/d  plot  would  be  met.  Essen¬ 
tially  our  recommendations  acknowledge  that  additional  experi¬ 
mental  data  could  be  beneficially  utilized..  We  are  attempting 
in  these  recommendations  to  direct  any  additional  experiments 
which  might  be  conducted  towards  specific  objectives  which,  when 
addressed,  will  provide  the  most  additional  insight  into  bomb  cra¬ 
tering  of  runways. 
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SECTION  VIII 


SUMMARY 

In  this  report,  a  model  analyaie  is  conducted  and  proce¬ 
dures  are  recommended  for  using  small  scale  experiments  to 
determine  the  extent  of  bomb  crater  damage  to  rigid  runways. 
Emplilcai  lorrr.uli  were  also  devised  for  estimating  true  crater 
size  and  the  extent  of  concrete  cracking.  Before  one  can  decide 
how  charge  weight  should  be  scaled  in  a  model  or  what  empirical 
formula  applies  for  predicting  true  crater  size,  one  must  deter¬ 
mine  if  the  mode  of  response  falls  in  the  camouflet  or  cratering 
mode.  Different  physical  phenomena  are  important  in  the  two  dif¬ 
ferent  response  modes.  Gravitational  effects  are  significant  in 
the  camouflet  mode,  but  not  in  the  cratering  mode.  Concrete 
pavement  strength  and  thickness  play  essentially  no  role  in  the 
cratering  mode  except  perhaps  for  radius  of  concrete  cracking; 
however,  they  influer.i  surface  phenomena  such  as  surface  crater 
radius  and  extent  of  cracking  in  the  camouflet  mode.  One  refers 
to  a  plot  of  W^*/d  versus  W^/d  where  W  is  charge  weight 

and  d  is  depth  of  burial  to  determine  whether  the  cratering  or 

1/4 

camouflet  mode  predominates.  If  the  value  of  W  /d  to  be  tested 

1/4  * 

is  greater  than  the  calculated  value  of  W  / d  from  Equation  (4), 

*  Equations  in  this  summary  are  numbered  as  in  the  text. 
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the  response  falls  in  the  camoufl'-t  mode,  and  if  this  value  is  less, 
it  falls  in  the  cratering  mode. 
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In  the  cratering  mode,  a  replica  modeling  law  is  appropriate, 
and  charge  weight  scales  as  the  geometric  scale  factor  cubed. 
Experimental  data  indicate  that  the  presence  of  the  pavement  has 
an  insignificant  influence  on  true  crater  size.  Empirical  equations 
obtained  by  curve -fitting  least  squares  straight  lines  to  nondimen- 
sional  pi  terms  from  a  similitude  analysis  yield  Equation  (5)  for 
predicting  true  surface  crater  radius,  Equation  (6)  for  predicting 
true  crater  depth,  Equation  (7)  for  predicting  true  soil  crater 
volume,  and  Equation  (8)  for  predicting  soil  j  us  concrete  crater 
volume.  One  standard  deviation  calculated  from 
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experimental  data  equals  essentially  12%  for  all  normalized  de¬ 
pendent  variables  except  R^/d  where  the  standard  deviation 
equals  approximately  20%.  Although  the  author  believes  this  ac¬ 
curacy  is  excellent,  ideas  are  suggested  which,  if  followed,  might 
improve  the  accuracy  moderately.  Radius  of  concrete  cracking 
in  the  cratering  mode  can  be  estimated  using  Equation  (9)  pro¬ 
vided  the  ultimate  strength  of  the  concrete  is  approximately 
1  0,  200  psi. 


=  4.247 
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The  amount  of  cracking  probably  depends  to  a  moderate  amount 
on  scaled  concrete  strength,  but  not  on  scaled  depths  of  burial. 
Because  most  of  the  data  are  on  concrete  pavements  with  strength 
of  10,  200  psi,  empirical  formulae  could  not  be  devised  for  other 
strergrh  concretes  in  the  cratering  modes.  Unless  concrete  strength 
cause  changes  in  mode  of  failure  or  energy  transfer  mechanisms, 
higher  strength  concretes  probably  crack  less  and  weaker  strength 
concretes  probably  crack  more  than  estimated  using  Equation  (9). 

In  the  camouflet  mode,  an  empirically  obtained  compro¬ 
mise  coefficient  is  used  to  scale  charge  weight  as  the  geometric 
scale  factor  to  the  7/24  power.  This  coefficient  is  a  compromise 
between  the  one  fourth  power  associated  with  scaling  gravitational 
and  inertial  effects  and  the  one  third  power  associated  with  scaling 
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strength  or  constitutive  and  inertial  effects.  Empirical  equations 
for  predicting  crater  size  in  the  camoufiet  mode  have  been  devel¬ 
oped  from  least-squares  fits  to  straight  lines  on  log-log  plots  of 
nondimensional  pi  terms.  Equation  (10)  is  for  predicting  maxi¬ 
mum  crater  radius,  Equation  (11)  is  for  predicting  true  crater 
depth,  and  Equation  12  is  for  predicting  true  crater  volume. 
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One  standard  deviation  was  predicted  for  these  curves  and  equaled 
approximately  12%.  Because  the  range  in  charge  weight  is  nowhere 
as  large  in  the  camoufiet  mode  as  in  the  cratering  mode,  this  low 
standard  deviation  is  not  as  significant.  Radius  of  concrete  cracking 
in  the  camoufiet  mode  can  also  be  estimated  provided  that  the 
strength  of  the  concrete  is  approximately  10,200  psi.  Equation  (13) 
is  the  equation  for  predicting  radius  of  cracking. 
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One  standard  deviation  for  R^/d  equals  11%  even  though  scaled 
thickness  of  concrete  and  scaled  slab  dimensions  are  not  consid¬ 


ered  in  this  analysis. 


This  report  closes  with  a  list  of  five  recommendations  for 
future  experiments.  We  suggest  that  additional  data  be  obtained 
to:  1)  establish  the  curve  separating  the  cratering  and  camouflet 
modes  over  a  wider  range  in  initial  conditions,  2)  measure  true 
crater  size  over  a  wider  range  of  charge  weights  in  the  camou¬ 
flet  mode,  3)  study  the  influence  of  concrete  strength  on  surface 
related  phenomena  in  both  response  modes,  4)  determine  the 
influence  of  soil  strength  especially  the  seismic  velocity  of  the 
soii  on  crater  size,  and  5)  evaluate  the  effects  on  bomb  detonation 
inside  the  pavement  on  crater  geometry  and  concrete  cracking. 

We  believe  much  additional  valuable  data  can  be  obtained  using 
model  experiments. 
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